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GRAVITATIONAL, AND DISSIPATIVE FORCES 
by FredS. Sidransky 
Lewis Research Center 

SUMMARY 

The method of characteristics is used to present general compatibility relations for 
nonsteady liquid flow or water-hammer theory which permit the investigation of the dy- 
namics of the flow under diverse conditions; that is, for a liquid flowing in a variable 
area conduit affected by heat addition, by gravity, by friction forces, and by centrifugal 
forces such as in a mixed- flow radial vane pump. The compatibility equations are pre- 
sented in finite difference form, and numerical-graphical methods for their solution are 
indicated. 


INTRODUCTION 

In rocket engine and Rankine cycle systems, the analysis of nonsteady liquid flow is 
often complicated by the transmission of the fluid in a conduit with a gradually varying 
cross-sectional area. In addition, the fluid may be affected by heat addition and a num- 
ber of forces such as gravity, dissipative, and centrifugal (e. g. , as in a pump). Hence, 
a general method is required which will allow the investigation of fluid transients under 
such diverse conditions. 

Analyses for nonsteady liquid flow have been presented in a number of papers for 
particular conditions; that is, for a tapered line (ref. 1), for friction (ref. 2), and for a 
|>ump as a lumped component in a conduit (ref. 3). 

In this report the theory of nonsteady liquid flow as developed in reference 4 will be 
extended so that the dynamics of a liquid in a conduit with a gradually varying cross- 
sectional area, undergoing heat addition, and affected by the previously mentioned forces 
may be analyzed. Applications will be outlined for nonsteady liquid flow with heat addi- 
tion, for tapered conduits with gravity and dissipative forces, and for the internal dy- 
namics of a mixed flow pump with radial vanes. 



ANALYSIS 


\ 


Derivation of Fundamental Compatibility Relations 


Continuity. - The continuity equation (ref. 5) is given by 


— (pA) + A (puA) = 0 (1) 

at as 


3 A 

and if — = 0, equation (1) may be expanded to 

at 

A&+(&\i+^p\ A+i^pu = 0 (2) 

at \as as / as 

3 A 

The term — , a function of the elasticity of the conduit, is omitted since it is usually 

g A 

negligibly small compared to the variation of area with distance along the conduit 

If equation (2) is divided by pA, then ^ 

lie+l u i£+i2»u 3 fo A >=o (3a) 

p at p as as as 


and since 


P=f(p, s) 


(3b) 


equation (3) may be expressed by 


i fifi. ap + i£ as\ + 1 /ag. a p + a s \ + — + U 5 ( ln A ) = 0 (4) 

p Vap at as at/ pVapas as as/ as as 


Assuming that 

as 

velocity, namely, 


is negligibly small for a liquid and using the definition for the acoustic 



( 5 ) 


result in 


2 


( 6 ) 
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Inasmuch as the pressure P is expressed in force per unit area and the density p in 
mass per unit volume and the head in the units of pressure divided by the specific weight 
of the fluid in the Earth's gravitational field, the following relation may be deduced 


P - P o 

P 


= Sc H 


(7) 


where P Q is the datum pressure, a constant. In view of equation (7) and the assumption 
that the variation of the density has a negligibly small effect on the head (i. e. , p = const. ), 
equation (6) may be presented by 


gc8H + g Cli aH t 8u + „ 8(lnA) _ 0 
a 2 8t a 2 8| 84 

2 

Multiplying by a /g and transposing some terms result in 

8H f T| 8H + a 2 8u _ ua 2 8 (In A) 

at 84 g c 34 ’ g c H 

Momentum . - Consider the momentum equation given by (ref. 5) 

iH + u iH = .ii p + f 

8t 84 p 8| 

Multiplying equation (10) by a/ g and considering equation (7) give 

in + i_ u in + a iS = f 
g c at g c 84 34 g c 

Subtracting equation (11) from equation (9) yields 

±/h-^u\ + (u- a) jg- -g-fu iS- *. f 

at V g c / a « g c V a «/ g c g c 


( 8 ) 


(9) 


( 10 ) 


( 11 ) 


( 12 ) 
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if a/g is a constant. By further algebraic manipulation, equation (12) can be pre- 
c 

sented in the form 

i./ H --2-u) + (u - a) — (h - — u) if 

9t y g c / V g c ' ®c ^ S C 


( 13 ) 


Adding equations (9) and (11) yields another expression, namely, 

JL (h + -g- u\ + (u + a) ± (h + £. u\ = - ^ 9(ln ^ f 
9t \ g „ / SI \ g c / g„ SI g c 


(14) 


Now let 


H + — u = C 


(15) 


and 


H - _ u = C_ 
Sn 3 


(16) 


and define two directional derivatives as follows: 


*! = i.+ ( u + a)i- 

6t at as 


(17) 


and 


£ = ± + (u-a)i- 


6t at 


SI 


(18) 


With the help of equations (15) to (18), equations (13) and (14) may be expressed by 


6 C y _ ua 2 3(ln A) a f 
St g c s| g c 


(19) 


5 + C_ 


6t 


\X3l d(]n A) + a £ 
Sr 9 I Sr 


(20) 
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^In the terminology of the method of characteristics (ref. 5, vol. 1, appendix A), equa- 
tions (19) and (20) can be considered as right running and left running compatibility equa- 
tions respectively, and their characteristic directions in the time-distance plane are 


*l = u- 
dt 


a 


( 21 ) 


along the right running characteristic and 


*i = u + a 
dt 


( 22 ) 


along the left r unnin g characteristic. The similarity of these results and those of refer- 
ence 4 is evident. 

If the flow velocity is very small compared to the acoustic velocity, equations (21) 
and (22) become simply 



dt 


(23) 


a 

dt 


(24) 


If the two terms appearing to the right of equations (19) and (20) are omitted or set equal 
to zero, then 


6 + C 


1 = 0 


fit 


(25a) 


and 



(25b) 


is the right running compatibility equation and 

C =H+— u 

e c 


( 26 ) 
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is the left running compatibility relation. These equations are the more familiar water 
hammer relations. 

Energy . - The energy equation gives rise to a third compatibility relation and char- 
acteristic direction. This may be seen by considering the energy equation in a conduit as 
expressed by (ref. 5) 


qpA as JL 

pa(s + A 

+ ± 

at 

L v 2 /j 



puAii + Pv + — 

V 2 ) 


(27) 


which after some manipulation becomes 


qpA = (i + 


u 


9(pA) , 3(puA) 

at as 


+ 9(hap) +pA fii + u9i) + pA 
a£ Vat as/ 


d_ 

at 


v 2 / as V 2 


By the definition of the co-moving or substantial derivative, namely, 

D^ = _a_ + ds _a_ 

Dt at dt as 


(28) 


(29a) 


or 


»=i- + u i_ 

Dt at as 

and, in view of equation (1), equation (28) may be presented by 

p MsA) + ^ clP + pA Dl + pAu Du = ^ 
as as Dt Dt 

From continuity, the following relation is obtained: 

d/‘) 
a(uA) _ Ap Vp/ 
as Dt 


(29b) 


(30) 


( 31 ) 


A A 

where — is omitted. Hence, instead of equation (30), 

at 


6 



J 

✓ 



is obtained. In view of equation (10), equation (32) may be transformed to 



Dt Dt 


q 


Since, from thermodynamics, it is known that 



it follows that 



( 32 ) 


(33) 


(34) 


(35) 


if the magnitude of the term uf , which represents in the negative sense the specific work 



Distance along conduit, £ 


Figure L - General characteristic net point for nonsteady 
liquid flow or water hammer. 


done by dissipative forces per unit time, is much 
smaller than the value of q (i. e. , the heat addi- 
tion is reversible). Equation (35) is the third 
compatibility equation; and in view of equations 
(29a) and (29b), the direction of the characteris- 
tic in the time-distance plane will be given by 

*Uu (36) 

dt 

Hence, in general, a net point in nonsteady liq- 
uid flow is defined by three characteristics, the 
projection of which on the time-distance plane 
is seen in figure 1. 

Heat addition. - Equations (15) and (16) are 
two equations with two unknowns, H and u, 
since C x is a constant along the left running 
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characteristic and is a constant along the right running characteristic, and the 
acoustic velocity a is also constant. If none of the parameters of equations (15) and (16) 
are affected by entropy change, which can be reasonably assumed in liquids, the contri- 
bution of heat addition to the dynamic analysis of liquids may be omitted. Yet it is desir- 
able to estimate the temperature variation of a liquid as it proceeds, for example, 
through a boiler or heated pump inlet so as to predict the incipience of boiling. This may 
be accomplished by following a particle of fixed identity in the fluid as implied in equa- 
tion (35). From thermodynamics, it is given that for constant volume reversible heating 

1 


T Ds = C y DT 


(37) 


Integrating equation (37) along the particle path (cf. fig. 1) results in 


T 3 - T 4 e 


(s 3 -s 4 )/C l 


and, consequently, by equation (35) 


T 3 = T 4 e 




(38) 


(39) 


where 


* = <L 

T 

Evaluation of the Specific Force 

The specific force f (or force per unit mass) appearing in the fundamental momen- 
tum equation (eq. (10)) is undefined. It may be equated to a specific gravity force, or a 
dissipative force, or a specific centrifugal force, or the sum of all these depending on 
the particular study to be made. In the following, the specific force f will be equated, 
first, to the sum of the specific gravity and specific dissipative forces and, then, in order 
to analyze the internal dynamics of a mixed-flow radial pump, to the sum of specific cen- 
trifugal and specific dissipative forces. 

Gravity and dissipative forces . - If the gravity force per unit mass, or in other 
words the specific gravity force g, is opposite in direction to the assumed positive direc- 
tion of the specific force f , the sum of the specific gravity and dissipative forces may be 
expressed by 
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( 40 ) 


f = - U U 
2D | u | 

where X is a friction coefficient and D is the hydraulic diameter of the duct. Hie fac- 
tor u/{u| ensures that the dissipative force is opposite in direction to the flow velocity. 
By equation (40) the compatibility equations (eqs. (19) and (20)) become 


6 *"y _ ua^ 3(in A) t a 

* g c H g c 

6 + C 2 a/ , A v 

x _ ua 3 (In A) a 

« g c 31 g c 

or, in finite difference form, 



/ Xu^ u 
\2D |«| 



< o 

Xu u 

2D l« 



\ Au^ u 

S c 2D M 


+ g 


I ( 4 3 “ 

23 


/ a Xu 2 u 

\g c 20 H 


+ g (*5 


13 


h* 


(41) 

(42) 

(43) 

(44) 


and, for the corresponding characteristic directions, 


*3 



V*2 


= “ a 23 


(45) 


*3 



a 13 


(46) 


where the subscripts refer to net points in the characteristic network (cf. fig. 1). The 
subscripts 13 and 23 represent average values between net point 1 and 3 and net 
point 2 and 3, respectively. Clearly, equations (43) and (44) may be more easily eval- 
uated if the variation of area with conduit length is expressed in the form 
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Figure 2. - Mixed-flow pump with radial vanes. Seeing that 
r - r f * f; sin o, it follows that f ■= u?E; sin^ a + iAj sin o. 


vanes may be expressed by (cf. fig. 2) 


A = e h ^) 


(47a) 


3(ln A) = MO (47b) 

h n 

where h(| ) is some easily differentiable 
function of conduit length. 

Pump dynamics . - The equation of 
motion for a mixed-flow pump with radial 


0u 0u ( 2 t . 2 2 . \ 

— + u lu> £ sin a + a? r. sin a) = 

0t 0| \ 1 / 


where u is the through-flow component of the flow velocity relative to the blade and o 

is the angle the through-flow component makes with the axis of the pump. The terms 
2 2 2 

at £ sin a + w r^ sin a represent the centrifugal force per unit mass. If, in addition, 
friction forces are significant, then 

2 

f = co 2 £ sin 2 cr + oi 2 r. sin a - — — — 

1 2D I u I 


and equations (41) and (42) may be modified to 


Sr Sr 


!i = .^M.iL2 $sin 2 ff + A S in<r-*LJL 


2D u 


: x ua 2 0(ln A) a / 2 t . 2 2 

— = + — I co £ sin a + co r. sin a 

t g r H g r \ 1 


Xu 2 u 
2D | u 


or in finite difference form to 


C V,= C V„- ~ ^ «3-‘2> 


^3 »* [g c a? 


— sin^u + co^r. sin a - — — — 

g c \ 2D |u 


£nrr «*-*■> 
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c = c 

x 3 x 2 


ua _ 3(ln A) 


L g c 


(to - t«) + 


13 


O 

a [ 2 v . 2 2 a Xu u 

— [ co £ sin a + a> r . sin — 

g c V 1 2D |u| 


(t 3 - tj) 


J 13 


(52) 


If cavitation occurs in the pump (the head is less than the vapor head), a vapor region 
will appear in the characteristic network. The vapor head and the C x characteristic 
will define the flow into the cavitation zone, and the C y characteristic will define the 
flow out. Hie time increment and past conditions will help to estimate the augmentation 
in the size of the cavitation zone. Further details on the analysis of the dynamics of cavi- 
tation regions by essentially the method of characteristics may be found in reference 6. 


CONCLUDING REMARKS 

General compatibility relations have been presented which will permit the analysis of 
the dynamics of a liquid under diverse conditions. In view of their complexity, a high- 
speed computer should be employed in their practical application to any particular 
problem. 

Although an electromagnetic pump has not been considered in this memorandum , it 
is clear that in place of the centrifugal force per unit mass term in the compatibility equa- 
tions appertaining to pump dynamics, an electromagnetic force per unit mass term may 
be employed. 

The characteristic methods of this report may be applied, for example, to an approx- 
imate estimate of the effect of the through-flow channel of a pump on flow dynamics. Fur- 
ther, the techniques of this report may be used to predict the dynamic formation and col- 
lapse of cavitation regions in a heated pump inlet line as well as the incipience of two- 
phase flow in a boiler under dynamic conditions. 

Lewis Research Center, 

National Aeronautics and Space A dminis tration, 

Cleveland, Ohio, June 8, 1966, 

120-27-04-27-22. 
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APPENDIX A 


SYMBOLS 


A 

area, sqft 

s 

specific entropy, ft-lb/ 

a 

acoustic velocity, ft/sec 


(slug)(°R) 

c v 

specific heat at constant vol- 

T 

temperature, °R 


ume, ft-lb/(slug)(°R) 

t 

time, sec 

C x 

parameter (see eq. (15)), ft 

u 

flow velocity, ft/sec 

C y 

parameter (see eq. (16)), ft 

V 

specific volume, cu ft/slug 

D 

hydraulic diameter, ft 

X 

friction coefficient 

f 

specific force or force per 

A 

i 

distance along conduit, ft 


unit mass, ft/sec 

n 

p 

density, slug/cu ft 

g 

specific gravity force, ft/sec 

a 

angle that the relative 

g C 

acceleration due to Earth's 

o 

gravitational field, ft/sec 


through- flow component of 
velocity makes with pump 

H 

head, ft 


axis 

i 

internal energy, ft-lb/slug 

* 

parameter, q/T 

•a 

P 

pressure, lb/sqft 

0) 

angular speed, sec" 

q 

heat added per unit mass per 

Subscripts: 



unit time, ft-lb/(slug)(sec) 

o 

datum 

r 

radial height (see fig. 2) 

1,2, 3,4 

refers to net point locations 

r i 

inlet radial height (see fig. 2) 


(see fig. 1) 
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